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I. INTRODUCTION 

The dynamics of glassy systems is by definition 
out-of-equilibrium over experimental time-scales. 
Classiness manifests itself through a great vari- 
ety of dynamical features such as aging, non- 
exponential relaxation of correlation functions or 
super-Arrhenius slowing down of the dynamics 
However the definition of glassiness remains 
an open problem because, in many situations, no 
static parameters were found to indicate whether 
or not a system is in a glassy state. There is indeed 
a general agreement about the fact that the glassi- 
ness of a system does not necessarily arise from an 
underlying static transition 

To shed light on this question several dynami- 
cal approaches have been developed Here we 
explore the suggestion made in [3[ that glassiness 
arises in a system when a coexistence between ac- 
tive and inactive regions of space-time takes place. 
The idea is that these dynamical heterogeneities 
are a defining feature of glassy systems. A method 
to inquire into the space-time character of a system 
is the application of the thermodynamic formalism 
of histories developed by Ruelle and coworkers Q ■ 
While the equilibrium statistical formalism stud- 
ies the fluctuations in the configuration space of 
the system, Ruelle's formalism focuses on the time 
realizations the system can follow in configuration 
space. A central parameter in this formalism is 
the activity K(t) of a history, which is the num- 
ber of changes of configurations between the initial 
time, set to zero, and time t. This parameter is a 
physical, time-extensive observable. 

For systems without quenched disorder, it was 
suggested that when for large times a system is 
found to have two well separated sets of histories 
- one where the activity is extensive in the system 
size, and the other where the activity is subex- 
tensive - there will be a coexistence of active and 
inactive regions of space-time in the system, sep- 
arated by sharp interfaces. Thus the system will 
exhibit slow dynamics. This was indeed success- 



fully shown in the case of Kinetically Constrained 
Models of glasses Q and for the non-equilibrium 
steady state of some Markov processes [6[ . 

Here we investigate the space-time properties of 
a family of disordered systems. Besides the fact 
that these exhibit experimental features akin to 
those of structural glasses, there is some intrinsic 
interest to focus on the effect of quenched disor- 
der. As will become clearer in the next sections, 
the central part of our study bears on the statistics 
of the activity K(t), which is defined for each real- 
ization of disorder. From the probability distribu- 
tion of K one can build a mathematical object -its 
large deviation function- which in many respects 
can be argued to play the role of a dynamical free 
energy. As is known from the statics of disordered 
systems, the presence of quenched disorder induces 
additional difficulties when it comes to averaging 
over the disorder degrees of freedom. We have had 
to deal with similar ones for the dynamical free 
energy we are after in this work. 

In this paper we focus on the directed Trap 
Model and the Random Energy Model. Although 
they are simpler to study they reproduce exper- 
imental features of more complex disordered and 
glassy systems like the super-activated behavior 
of the viscosity in glasses Q , aging 0, Q or non- 
trivial violations of the fluctuation-dissipation the- 
orem In both cases, we show analytically 
that the aforementioned large deviation function 
displays a discontinuity between the sets of active 
and inactive histories : this suggests that also in 
these disordered systems there is a signature of the 
glassy behavior using the thermodynamic formal- 
ism of histories. Furthermore we study the influ- 
ence of the distribution of the disorder on the pres- 
ence of a dynamic transition and discuss its occur- 
rence in connection (ill . [l2| with the distributions 
of extrema of disordered energy landscapes. 
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II. DIRECTED TRAP MODEL 

A trap model is defined by N independent traps 
labeled by an integer i, each trap being character- 
ized by an energy Ei. We consider a continuous 
time Markov dynamics among traps: the dynam- 
ical evolution is specified by the probability P%{t) 
that the system stays in trap i at time t, and by 
the transition rates Wij for jumping from trap j 
to trap i. In the directed trap model the transition 
rates take the form 
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where 



B j = e 
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It can be shown that Pi(s,t) obeys an evolution 
equation (similar to but different from a Master 
Equation) of the form d t P l {s,t) = W K P l {s,t) 
where the elements of the evolution operator take 
the form: 
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From p(s,t) in Eq. @. one defines the cumu- 
lant generating function of the activity K: 
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i = 1,...,N. Thus the system evolves through the 
following Master Equation: 



dPi(t) P{t) P;_ x (t) 
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The trapping times Bj are random variables dis- 
tributed according to a Levy distribution 



p(Bj) = 9(B 3 - l)/xBT 



H G (0, 1) , (3) 



where 6{x) is the Heaviside step function. In 
the following we will use the notation O(Bj) = 
J p(Bj)0(Bj)dBj to denote an average over Bj. In 
particular, the distribution in Eq. ([3]) implies that 
the mean trapping time is infinite, which causes 
anomalous diffusion. 



A. Large Deviation Function 

Ruelle's thermodynamic formalism enables to 
investigate the dynamics of a system by studying 
the histories the system can follow in the configu- 
ration space. In the directed Trap Model configu- 
rations are represented by traps so that a history 
in the configuration space between time and time 
t is specified by a sequence of traps visited during 
this time and by the time intervals elapsed between 
each jump from one trap to the next. To per- 
form a statistics over histories one has to classify 
them according to a time-extensive (and history- 
dependent) parameter. A suitable one is the activ- 
ity K(t), defined as the number of configurations, 
here the number of traps, visited between time zero 
and time t. In the directed trap model, if the dy- 
namics starts at time t = in the trap i = 0, K{t) 
is simply the trap in which the system stands at 
time t. 

We now consider the probability Pi(K,t) to be 
in trap i at time t at fixed activity, and define its 
Laplace transform: 



where (...) stands for an average over all possible 
histories the system can follow between time and 
time t. At large time t one expects 



Z K (s,t) ~ e 



(7) 



where iPk(s) is the largest eigenvalue of the linear 
operator Wk m Eq. (|5| and is thus a large devi- 
ation function. At large time t, the derivatives of 
V'if ( s ) wm §i ve the cumulants of the activity K(t). 
Note that ipx{s) depends a priori on the realiza- 
tion of the random variables Bi §3§ and will thus 
be itself a random variable. 

The parameter s allows to probe the different 
histories the system can follow: positive s will 
correspond to the inactive phase of the dynam- 
ics, i.e. K smaller than its average, while negative 
s correspond to the active phase of the dynamics, 
i.e. K larger than its average. A discontinuity in 
the derivatives of the large deviation function will 
correspond to a dynamic transition between two 
different phases (the active and the inactive one), 
reflecting the sharp interfaces between dynamical 
heterogeneities. 



B. Dynamic Transition 

To compute the large deviation function i^k{s) 
we have to find the largest eigenvalue of W K 
defined in Eq. ©. The N eigenvalues Xj(s), 
1 < j < N ar c solutions of the equation 

_ s Pi-x(s,t) P(s,t) - 
e — — = \j{s)Pi(s,t) , (8) 



B, 



and ipKis) = maxi<j<jv {Aj(s)}. The computa- 
tion of the characteristic polynomial is straightfor- 
ward and leads to the eigenvalue equation 



1 N 

-5>(l + Ai(s)Bi) = -* 



(9) 



Pi(s,t) 



E 

K=0 



e- sK P(K,t). 



(4) Obviously, one has ^jk(s = 0) = 0. Let us first 
consider the case s < where there is only one 
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Xj > solution of Eq. © which thus coin- 
cides with ipx(s)- In that case, i.e. in the active 
phase, it is natural to assume that ^k{s) is self- 
averaging so that, in the limit N — > oo, one expects 
iPk{s) — ^Pk{s). Using the law of large numbers 
to treat the sum in Eq. ([9|), one obtains in the 
N — ► oo limit 
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dBB- 1 -»\n(l+~4^{sjB) 
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which uniquely determines iPk(s). To describe the 
fluctuations of i/>k(s) around its mean for finite 
N > 1 one writes i/) K (s)-ip K (s) = Xk(s)N~ 1 / 2 + 
0(N~ 1 ) where Xk{s) is a zero-mean random vari- 
able. Using the Central Limit Theorem applied 
to the sum in Eq. (|9|), one obtains that xk(s) is 
a Gaussian variable such that the distribution of 
iPk{s) is given for N ^> 1 by 



p(^k{s)) = — ^=exp 



where er^, is given by 
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x[ln"(l + 4>k{s)B) - (ln(l + 4>k{s)B)Y 



where ^k{s) is given in Eq. (|10p . 

Let us now consider s > and label the B^s 
such that < B\ < ■ ■ ■ < B^. In that case, one 
obtains the bounds 



— < if> K (a) < . 
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From the distribution of the variables BiS in 
Eq. one obtains that the one of the largest one 
.Bat is given by: 

p(B N ) = N»BJ?-»(1 - B^f- 1 , (14) 

from which one gets for N ^> 1: 
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Thus, from the bounds in Eq. (|13|) we find 
that liniTv^oo tJjk(s) = for s > 0. On the 
other hand, for s < 0, one has from Eq. (fTU|) 
liniA/^oo -0^ (s) > and in that limit iPk(s) oc 
(— s) 1 /' 1 for small s. Therefore, the first derivative 
tp' K {s) is continuous but higher order derivatives 

tpj(\s) with p > 1/fx will display a discontinuity 
in s = 0, indicating a dynamic phase transition (of 
order higher than one). 

We have checked the presence of this dynamic 
transition by solving numerically the eigenvalue 



equation ©. In Fig. [TJ we show a plot of 4>k (s) as 
a function of s for different values of N = 20, 30, 50 
and 100 for /j, = 0.8. The data were obtained by av- 
eraging over 10 6 samples. For s > these numeri- 
cal data show that ipn(s) — > when ./V — > oo (one 
finds indeed iPk(s) oc —N^ 1 ^, consistent with the 
bounds in Eq. (fl~3|) ). On the other hand, this plot 



on Fig. [T]shows that lim^v^oo iPk(s) > for s < 0. 
The so lid line in Fig. Q] is the analytical value of 
">Pk( s ) obtained by solving numerically Eq. (JTDJ) , 
which is in good agreement with the numerical 
data and shows the presence of a dynamical phase 
transition. 

The slope of the large deviation iPk(s) function 
is related by Eqs. (|6|7p to the mean value of the 
activity (K) over all histories, thus in the active 
phase (s < 0) of the dynamics the activity takes 
a constant value (A") > while in the inactive 
phase (s > 0) it is found to be subcxtensive in 
the system size ((A) ~ N^ 1 ^). This coexistence 
of active and inactive phases of space-time, or dy- 
namical heterogeneity is argued to manifest itself 
through the glassy properties of the directed Trap 
Model. In the corresponding pure system, where 
the waiting times are no longer random variables, 
glassincss obviously disappears, and so does the 
discontinuity in the derivatives of the large devia- 
tion function. The activity is a constant (A') > 
also in the s > phase and no dynamical hetero- 
geneity is present. 
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Figure 1: Plot of ipfc(s) as a function of s for different 
values of N and /j, = 0.8. The solid line is the solu- 
tion of Eq. flQj: for s > the curve approaches the 
horizontal axis as N' l/ ^ for N > 1. 

The numerical results in Fig. [1] were obtained 
assuming that the large deviation function ap- 
proaches its average in the large N limit and 
is distributed according to Eqs. (|11I12|1 . In- 
deed this supposition was confirmed evaluating nu- 
merically the probability distribution p(ipx(s)) of 
"ipK(s) for a negative s. In Fig. [2] we show a plot 
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of p{ip K {s))N~ 1 / 2 as a function of N l l 2 ($ K (s) - 
ip K {s)) for p = 0.8 and s = -0.4 (N = 100,200). 
The data were obtained by averaging over 3 ■ 10 7 
samples. The solid line is the Gaussian distribu- 
tion obtained rescaling by s/~N the Eqs. (fTTl [T2|) . 
This gaussian id approached by the numerical data 
as the system size grows. 

On the other hand, for s > 0, the evaluation of 
the average large deviation function, i.e. the solu- 
tion of Eq. (JTUJ) , was done limiting the values of 
tpK(s) to the iV-dependent bound (flU]) . Indeed we 
found numerically that the distribution p{ipK{s)) 
of iPk(s) for a positive s is well fitted by the form 
(US: 



p(ip 



>k(s)) = " ^k{sY) n ~' 



(16) 



for all N. In the inset of Fig. [2] we show the 
rescaled distribution p{i[>k{s))N~ 1 i function 
of ^(s)TV 1 ^ for N = 100, 200 averaged over 10 7 
samples for s = 0.2 and p = 0.8. These distribu- 
tion are found to be well approached by the large 
N limit of the rescaled distribution which takes the 
form p(x) = px^^ 1 cxp (— x^). 
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Figure 2: Numerical evaluation of the distribution of 
the large deviation function for p = 0.8 and s = —0.4 
(N=100,200). The solid line is the rescaling of the 
Gaussian distribution (JTTJ . Inset: numerical evalua- 
tion of the distribution of the large deviation function 
for p = 0.8 and s — 0.2. The solid line is the large N 
limit of the rescaling of distribution (|16[) . 



A. Dynamics of the REM 

The Random Energy Model is a system of 2^ 
energy levels Ei which are independent random 
variables distributed according to a Gaussian dis- 
tribution of zero mean and variance ^AfJ 2 : 



P(E) 



e -E 2 /ATJ 2 



(17) 



The model exhibits a phase transition at the crit- 
ical temperature T c = J/2v4n 2 below which the 
free energy becomes a constant Eo(J\f) ~ MJ\/hi 2 
in the thermodynamic limit and the entropy van- 
ishes. As a consequence, the energy levels close to 
Eq(M) will dominate the low temperature phase. 
We thus consider a system of N energy levels 



Ei =£b(A0 + < 



(18) 



where is small and nonextensive [l4| distributed 
according to 



0, 



e > e c 
e < e c 



where p = TjT 1 is a constant and e c is a cut-off 
energy. In the following, we will consider the scal- 
ing limit N — > oo, e c — > oo, while Ne~ pec = v 
is kept fixed, where the physical properties of the 
REM for temperatures T < T c take finite, cut-off 
independent values. 

The dynamics of this model is defined by a Mas- 
ter Equation 



dPjjt) 
dt 



^uv ; ;/i Y^\\„r,:n. (19) 



where the transition rates Wji for going from level 
i to level j are taken of the form [lj] : 



u„ = /;,r,r, 



with Bj = , (20) 



where i = 1 , . . . , N and Vi arc positive random 
variables. One can easily check that the rates Wji 
obey the detailed balance condition. The distri- 
bution of the BiS can be easily derived from the 
distribution of the energies e,;, leading to: 



for (#)""< B, 
otherwise 



(21) 



III. RANDOM ENERGY MODEL (REM) 

We consider now a widely explored model of dis- 
ordered s yste m, the Random Energy Model intro- 
duced in [l3| as a mean field spin-glass model. Its 
dynamics has been studied in detail by Kopcr and 
Hilhorst in Ref. El. 



where we set u = p/(3 = T/T c , p £ (0,1). Fol- 
lowing Ref. [14j |. we will consider level dependent 
barriers Vi = B~ q with q e (0, 1). If one associates 
to each level i an independent random magneti- 
zation pi with zero mean and variance N, the ki- 
netic REM can be considered as a magnetic model. 
Within this choice of level dependent barriers it has 
been shown that the behavior of the equilibrium 
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autocorrelation function for the magnetization ex- 
hibits a non-exponential behavior [141 ] 

C(T) cx r 2 ^ /q for q + pi > 1 , (22) 
C(t) oc t-^e- Dtl for g + m < 1 and g < - , 

C(t) cx £-W(2</-i) for g + /u < 1 and g > ^ , 

where D is a constant, 7 = fi/(l — 2q + fiq/ (1 — g)) 
and 77 = (2/i/(l - g) - 1) 7 . 

Here we apply the thermodynamic formalism of 
histories to this model and show that a dynamic 
transition takes place reflecting the slow dynamics 
in this system (|22]) . 



B. Dynamic transition 

To investigate the presence of a dynamic tran- 
sition, we look again at the large deviation func- 
tion i(jk(s) defined above in Eq. ([7]), associated to 
the number of configuration changes K(t) between 
time and time t. From the Master Equation gov- 
erning the dynamics (Tl9|) . one obtains the equation 
of evolution of Pj(s, t) defined in Eq. Q 



(23) 



= e- 3 BiVi VjPj(s, t) - ViPi(s, t) BjVj . 

The large deviation function ipj((s) is the largest 
eigenvalue of the evolution operator W . The 
eigenvalues Xj (s) of W K are solutions of the equa- 
tion 



/(A) = E 



e~'BiV? 



(Vi + (e- s - l)B t V: 2 + A 



(24) 
(25) 



where iPk{s) = rnaxi<j<jv {Aj(s)} and £ = 

We first focus on the case q+fi > 1. One notices 
that /(A) has TV simple poles on the negative real 
axis at x t = -((Vi + (e~ s - l)BiV 2 ) < 0, and de- 
creases to zero for A — ► 00. One has also /(0) = 1 
for s = 0, whereas /(0) > 1 for s > and /(0) < 1 
if s < 0. Thus in the active phase, s < 0, there 
is only one positive eigenvalue Xj which coincides 
with %Pk(s) > 0. As done previously for the di- 
rected trap model, we thus suppose that i^k(s) is 
self averaging in the limit N — ► 00 and coincides 
with its average. One obtains that, for s < 



il> K {s) ~ N^+f-V^ipKis) 



where iPk( s ) is independent of N and solution of 
the equation 



fxv l = dB = e s 

v Vt> vB-i + ipx(s) 



(27) 



where v = [w( 1 ~ q '' 11 /(q + fi — 1) > 0. Performing 
an analysis similar to the one done for the directed 
trap model, one also finds th at the fl uctuations of 
iPk{s) around its mean value iPk(s) are described 
by a Gaussian as in Eq. (fTTj) of variance iV" 1 / 2 . 

However, the analysis in the inactive phase is 
more subtle. Indeed, given the singularities of /(A) 
in Xi, together with the fact that /(0) > 1 for 
s > 0, iPk(s) satisfies the bounds 

- min(C^ + (e- s - l)B t V 2 ) < i' K (s) < . (28) 

i 

If one labels the random variables Bi such that 
B\ < B2 < ... < Bn one easily sees that the mini- 
mum in the left hand side of the inequality (f28|) is 
reached for the maximum Bn, whose distribution 
is given in the large TV limit by 



p(B N ) = v/iB 



N 



y-n e -vBi 



(29) 



Using this distribution (j2"5|) together with the dis- 
tribution of the Bi 's in Eq. (f2Tj) one obtains that 



ci 



N ( q+ll -i)/^ < q j jK ^ < Q ^ ^ 



where c\ > is a constant, independent of N. 
To find an upper bound for ipj((s), one writes the 
equation above (|24|) in a different form: 



E 



e- s B t V t 



e~ s C + (1 - e- s )(C - BiVi + A(s)Ur 1 ) 



= 1 . 



Since the first term in the denominator is the 
sum of the numerators, the remaining terms (1 — 
e~ s )(C — BiVi + X(s)V~ 1 ) must be negative for at 
least one i. In this way we get an upper bound for 

%l) K {s): 

lM«) < (e- s - 1) min(CK i - B t V 2 ) (31) 

i 

Again, the minimum in the right hand side of this 
inequality is reached for the maximum B jv. Per- 
forming the average, one gets 



(32) 



where < ci < c\ is a constant. Combining these 
two bounds (f30| l32j) one gets for s > 

-iMs) =0{N^ + ^^) . (33) 



(26) 



Thus by looking at the behavior of i/jk{s) for large 
N in Eq. ((251 S3 , onc shows that there is a 
dynamical transition as s crosses 0. Eq. ([2"7]) shows 
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that 4>k{ s ) <x ( — s) 9 /( 9 +^ -1 ) for small s indicating 
that, as in the directed trap model, the order of 
this dynamical transition is larger than one. 

These analytical predictions for q + ji > 1 in Eq. 
(|261 I33| have been confirmed by solving numeri- 
cally the eigenvalue equation (|24|) . The average 
value tPk(s) was computed by averaging over at 
least 10 6 realizations of the disorde r. In Fig. [31 we 
show a plot of N~^ +2q ~ 1 ^ ^i'K (s) as a function of 
s. In agreement with Eq. (|26[ I31?]) . one observes a 
dynamical transition occurring at s = 0. Moreover 
it has been checked that for positive s the large de- 
viation function behaves like predicted in Eq. (1331) , 
i.e. -i}jk{s) = 0(AT(9+m-i)/m). 

We have checked numerically that this transi- 
tion is also present for q + /i < 1, where the 
equilibrium correlation function displays a non- 
exponential decay (f22j) . In Fig. [U one shows a 
plot of N~ 1 / 2 ipK{s) as a function of s for <j = 0.45 
and (j, = 0.45. In that regime, the autocorrela- 
tion function decays like a stretched exponential 
and our numerical data clearly shows a transition 
occurring for s = 0. Similarly, in Fig. [5l one 
shows a plot of %l) K (s)N^ 2 ^ q ~ tl ^ tl as a function of 
s for q = 0.6 and = 0.35, where the correlation 
function decays like a power law. Here again, our 
numerical data clearly shows a dynamical transi- 
tion occurring for s = 0. Again these results were 
obtained averaging over at least 10 6 realizations 
of the disorder. Thus in the three cases we have 
an anomalous behavior of the correlation function 
(f2"2"j) and a discontinuity in the derivatives of the 
large deviation function i/jk(s) at s = (Fig. [3J[5]). 
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Figure 4: Numerical evaluation 
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Figure 3: Numerical evaluation of the large deviation 
function tpK(s) for q = 0.9 and fi = 0.8 (q + fj, > 1) in 
the Random Energy Model. 



Figure 5: Numerical evaluation of the large deviation 
function for q = 0.6 and /i = 0.35 (q + fj, < 1 and 
q > 1/2) in the Random Energy Model. 



IV. CONNECTION WITH THE EXTREME 
VALUE STATISTICS 

The glassy behavior of disordered systems 
emerges in the low temperature phase or in the 
long time limit. Consequently the extreme values 
of the disorder realizations play a dominating role 
with respect to the typical values. The low tem- 
perature phase of the Random Energy Model is 
indeed described through the distribution of the 
lowest energy levels ((TBj) . which is found to belong 
to the Gumbel universality class of the minimum 
of variables which are unbounded and have a dis- 
tribution that decays faster than any power law to 
— oo. 

It was pointed out in [ll|, EH that the Gumbel 
class corresponds exactly to the one step replica 
symmetry breaking solution in the replica lan- 
guage. It was also highlighted that for the Weibull 



class, i.e. the extreme value statistics of bounded 
variables, and the Frechet class for power-law de- 
caying variables, the replica method cannot be 
used without substantial modifications. 

In order to shed light on how the extreme value 
statistics of disorder affects the glassy properties of 
the system in the low temperature phase we con- 
sider two versions of the Random Energy Model. 
We choose for the equilibrium distribution of con- 
figurations in the low temperature phase, i.e. the 
Boltzmann factors (f^Oj) Bj = e - ^ , both an ex- 
ponential and a uniform distribution. We find a 
dynamical transition in the first case but not in 
the second. These results are corroborated by a 
study of the equilibrium autocorrelation function 
in both cases. 



Exponentially distributed Boltzmann 
factors 






N=500 + 
N=200 x 
N=100 * 
N=10 □ 




***** 

*** 

** 



-0.04 0.04 0.C 



Figure 6: Numeri cal eva luation of the average large 
deviation function "4>k{s) for exponentially distributed 
Boltzmann fa ctors, w ith p — q — 0.9. Inset: numerical 
evaluation of -^k(s) for uniformly distributed Boltz- 
mann factors. 



We recall that the dynamics of the Random En- 
ergy Model can be described through the following 
Master equation: 



dPjjt) 
dt 



where the transition rates for level-dependent bar- 
riers are Wj% = N~ 2 PB^ q Br q : with p > and 

? 6(0,1). 

We define as a version of the Random Energy 
Model a system whose dynamics is still described 
by Eq. (|34|) . with the difference that the statistics 
of the energy barriers in the low temperature phase 
has changed. Instead of a Levy distribution we 
consider a probability distribution function of the 
Boltzmann factors of the form: 



p{B i )=p{e-< 3ei )=N6(B i )e 



— NBi 



(35) 



where the scaling in N has been chosen a posteriori 
to find a well defined thermodynamic limit of the 
equilibrium autocorrelation function (see section 

uTTTT]) . 

The large deviation function iPk(s) is again the 
largest solution of the eigenvalue equation 



E, 



S B 



l-2q 



C^ + (e-«-l)^- 2 " + Vic( S ) 



(36) 



We solved numerically the equation ()36p and aver- 
aged the solution over I0 7 realizations of the dis- 
order. The results are plotted in fig. [51 

As in the directed Trap Model and the Random 
Energy Model we observe a dynamical transition. 
For s < the large deviation function appears to 
be constant in the system size. For s > the curve 
approaches the horizontal axis for N ^> 1. 



Thus even if the moments of the Boltzmann fac- 
tors (or the mean trapping times) arc well defined, 
the system presents dynamical heterogeneities and 
we expect an anomalous slowing down of the dy- 
namics. 



B. Uniformly distributed waiting times 

Let's consider uniformly distributed Boltzmann 
factors: 



p{Bi) = Ne{B, 



N 



a + 1 

N 



- Bi , (37) 



where a is a finite constant with ^ [a, a + 1]. If 
we allow the B^s take the value B = 0, after some 
time the system will be frozen in a single energy 
level. 

Evaluatig numerically the large deviation func- 
tion of this system we find a complete indepen- 
dence of the large deviation function from the sys- 
tem size for all s (see the inset in fig. [6]), and no 
discontinuity of the large deviation function. Since 
we do not observe any dynamic transition in the 
space of histories, we expect this system to be dy- 
namically homogeneous. 



C. Correlation function 

To compute the equilibrium autocorrelation 
function of the Random Energy Model with ex- 
ponentially and uniformly distributed Boltzmann 
factors we recast the Kopcr and Hilhorst's study 
[HI of the relaxation properties of the Random En- 
ergy Model. 
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Starting from the Master Equation with 
transition rates (|20p . the authors proved that the 
autocorrelation function C(t) of any physical quan- 
tity that does not couple to the energies of the 
REM can be written as a complex integral of the 
form: 



C{t) 



ds e 



-stc 



N 



V 3 B 3 V k B k 



k^J {s-V 3 0(s~V k O 



2iri 



(38) 

with Z — X)i=i Bi an( i where the line integral is 
taken over a contour oriented counterclockwise and 
encloses all the poles of the integrand on the pos- 
itive real axis. If we average this expression, as in 
[lij], over a Levy distribution we find the asymp- 
totic behaviors (|22[) which characterize the glassi- 
ness of the Random Energy Model. 

Here, to corroborate the presence or absence of 
the dynamical transition obtained previously, we 
average Eq. (|38[) over an exponential distribution 
and then over a uniform distribution. 



where a is a constant. After some manipulations 
left in the Appendix lAl one has 



C(t) 



a+l 



d.r 



a+l 



, xy -t 



(44) 



The minimum of the integrand is reached for x = 
y = a and the maximum for x = y = a + 1, thus 
the integral is bounded from above and below by 
an exponential and we conclude that: 



C(t) 



-tc(a) 



(45) 



where c(a) is a positive function of a. Eq. (|45| is 
the asymptotic behavior of the correlation function 
of a non-glassy system. 



CONCLUSION 



1. Exponential distribution 

We consider exponentially distributed Boltz- 
mann factors 



p(B) = N0(B)e~ NB 
and transition rates of the form 



Wi.. 



N~ 2q B}~ 9 B7 9 



(39) 



(40) 



where the N dependence is chosen to give a well de- 
fined thermodynamic limit of the correlation func- 
tion. After some manipulations which details are 
given in Appendix [5] one finds 



C(t) = t^+T 



dx I dyxye 
Jo 



-i5TTg(x,i/) 



, (41) 



where g{x,y) = x + y + T(2 - q)/{x q + y q ). One 
easily finds that g(x,y) has a single minimum in 
x K+ for x* = y* = (<?r(2 - q)/4) 1 /^+ 1 ) and 



d 2 x g{x\y* 



d;g(x*,y* 



Thus the large 



time behavior C(t) is given by a saddle point cal- 
culation yielding 



C(t) ~ 27rtiTT x *V t5+1 tt*"*"> 



for t > 1 (42) 



which shows, given that q > 0, that the relaxation 
is indeed slower than exponential. 



The central statement of our work is that when a 
system presents glassy features (anomalous slowing 
down of the dynamics, non exponential decay of 
the correlation functions, anomalous dependence 
of viscosity on temperature and so on) it appears 
heterogeneous in time and energy landscape, and 
a dynamical transition in the space of time realiza- 
tions the system can follow in configuration space 
has to be found. 

Here we validated this statement for two well 
known disordered models of glasses, the Random 
Energy Model and the Directed Trap Model, de- 
tecting in both a dynamical phase transition. In 
the Random Energy Model the influence of the 
statistics of disorder over the space-time hetero- 
geneity of the system was analyzed. Indeed we 
found a dynamical transition for a REM with 
exponentially distributed Boltzmann factors and 
showed that we are in presence of a glassy dynam- 
ics by computing the equilibrium autocorrelation 
function of this system which is shown to decay 
anomalously, i.e. to be a stretched exponential in 
time. Furthermore we considered a uniform dis- 
tribution of Boltzmann factors finding no dynam- 
ical transition and showing, as expected, that the 
equilibrium relaxation is exponential in time and 
no slowing down of the dynamics is present. 



2. Uniform distribution 
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Appendix A: COMPUTATION OF THE 
CORRELATION FUNCTIONS FOR THE 
REM 



The equilibrium autocorrelation function ex- 
pressed in terms of these Green functions is 



In this appendix, we present some details con- 
cerning the computation of the corrrelation func- 
tion for the kinetic REM. 



a. General framework 



N 



,Ct( 



(A5) 



We recall the Master Equation (|34|) with transi- 
tion rates (PU1) : 



N 



dt ^ 

0=1 



VjPjM-ViPiMC. (Al) 



The probability to find the system at level i at time 
t if it was in j at time t = 0, i.e. the Green func- 
tion, can be written in terms of the eigenvalues and 
eigenfunctions of the Master Equation. We call <ft 
the ith component of the right-hand cigenfunction 
with eigenvalue —A: 



Using IjAip one can write: 



N 



ViBi Yuri 



(A2) 



(A3) 



where P cq = lim^oo Pj(t) oc Bj. 

To get rid of the dependence of the autocorre- 
lation function C(t) in Eq. (|A5j) on the unknown 
eigenvalues A; we consider its Laplace transform 
C(s), which using (|A4[) can be written as: 



C( S ) 



1 N 



if(A,) 



with Z = Bi and where 



Zf-^ ( S + A ; )F'(A Z ) 



(A6) 



N 



V f B o 



F(z) 



N 



The detailed balance condition implies that the 
left-hand eigenfunction is ip£ = (fif/Bi, so that the 
Green function takes the form: 



N 



G«(*)=£e" 

N 



Ait 



1=1 



VjBj Vj 



(A4) 



We introduce = H(z)/((s + z)F(z)), which 
has poles in z = — s, in z = Vj£ and in the zeroes 
of F(z). Looking at the eigenvalue equation (|A3[) 
we see that the sum of the residues of g(z) in the 
zeroes of F(z) will be exactly C(s). Thus, since 
the contour integral of g(z) along a circle centered 
in the origin with radius R vanishes as R —> oo, 
one finds: 



1 I E^LiE 



C(s) 



V J B J V k B k 
k ^0 (V j C+s)(V k C+s) 



V,C+s 



1 

Z 2 



N 

E^ 2 



(A7) 



o=i 



from which we see that C(s) has poles only on the 
negative real axis (s = is not a pole). Applying 
an Inverse Laplace transform and changing s into 
— s£ we obtain the expression of the autocorrela- 
tion function C(t) given in the text in Eq. (|55|) . 



b. Exponential distribution 

We consider exponentially distributed Boltz- 
mann factors 

p{B) = N9{B)e~ NB , (A8) 
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Using for the denominator in Eq. (f58|) the integral 
representation a" 1 — f Q d\e~ aX , we obtain from 



Eq. (|3o|) the expression for the averaged autocor- 
relation function: 



C(t) = lim A 4 



OO 



OO 



ds 



BiBj 



' ' dX J d »f 2vis(sB!-l)(sB]-iy 



-N(Bi+Bj) 



(A9) 



-X^Bi+Bj)-^- 



x e 



-stN- 2q {B]- q +B)- q ) 



X J 



iV-2 
N 



where the integral In is given by 



In = N I dBe-NB-XB-n^-rtN-XB*-' 



Thus using this large N expansion (j A10|) in the 
formula (|A9|) . and integrating over A and \l we get: 



Rescaling the variable B = B'/N and integrating 
by parts we have, to leading order in N: 

Jjv-I-— / dBe- B (X-^ + stN- q (l-q)B- q ) 
A Jo 

(A10) 



/v 



C(t) = lim N / dBi / dSj 



ds 



B{Bj 



-N{B i +B j )-stN- 2q (B 1 i - q +B 1 - q +N q V(2-q)) 



2ni s(Bi + Bj + 1) s(B q B 3 + B t B q ) - B t - B 3 - {sB q - \){sB q - 1) 



Changing variables according to NBi — x, NBj = text (l4Tj) . 
y and keeping only the leading orders in A we have: 



poo poo p j 



IQ JO 

e -(x+v)-atN-T(2-q) 

s(sN-*(x* + yi) - 1) ' 



(All) 



which has one pole in s = N q / (x q + y q ). Complex 
integration thus gives: 



C(t) 



dx 



dyxye-^+^e'*^^ (A12) 



c. Uniform distribution 

We now consider uniformly distributed Boltz- 
mann factors: 



p(B)=N6(B-^)9(^±-B), (A13) 



Finally we change variables x = y = 

y't^/ii+^j anc ] g nc i t ne expression given in the 



where a is a constant. The average of the correla- 
tion function (l38l) takes the form: 



J 



C(t) = lim N 

N^oo 



q+1 

N 



dB q 



N 



q+1 
N 



dBi 



N 



ds 



BiBj 



A /o **! 2ris(sB q -l)(sB q -l) 



x e 



-X(B t + B J )-^ 7 ^- T + 7 ^ j L- T )-stN- 2 "(Bl- q +B 1 ] - q ) 



-N-2 
N 



r 



(A14) 



where 



q+1 

Jn= l N dBNe- XB -^^r- stN ^ PBl - q 
"'77 



Since the interval of integration is very small when 
iV > 1 we can consider the integrand as a constant 
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over this interval and make the approximation: 



Jn 



-\aN- i + t j,aN--stN 



- 2p -(i- q ) l-g 



(A16) 



The integrand has one pole in s = N q / (x q + y q ) so 
after complex integration and choosing p = q we 
obtain the formula given in the text 



Thus after integration over A and \i and changes of 
variables x = NBi and y = NBj, the correlation 
function becomes: 



C(t) = lim 



a+1 

dx 

tN -2 V + q a l- q 



a+1 



dy 4> ir-.xy 

2m 



sa 2 (sN~ q (x q + yi) - 1) 



(A17) 
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